SPECTRAL STABILITY FOR SUBSONIC TRAVELING PULSES OF THE 

BOUSSINESQ 'ABC SYSTEM 



SEVDZHAN HAKKAEV, MILENA STANISLAVOVA, AND ATANAS STEFANOV 



Abstract. We consider the spectral stability of certain traveling wave solutions for the Boussi- 
nesq 'abc' system. More precisely, we consider the explicit sech 2 (x) like solutions of the form 
{ip{x — vrt),ip(x — wt) = (ip, const. ip), exhibited by M. Chen, [7], [8] and we provide a complete 
rigorous characterization of the spectral stability in all cases for which a = c < 0, b > 0. 



1. Introduction and results 



1.1. The general Boussinesq 'abed' model. In this work, we are concerned with the Boussi- 
nesq system 



The first formal derivation for this system has appeared in the work of Bona-Chen-Saut, [5] to 
describe the (essentially two dimensional) motion of small-amplitude long waves on the surface 
of an ideal fluid under the force of gravity. Here, r\ represents the vertical deviation of the free 
surface from its rest position, while u is the horizontal velocity at time t. In the case of zero 
surface tension r = 0, the constants a, b, c, d must satisfy in addition the consistency conditions 
a + b = t;(9 2 — 1/3) and c + d = ^(1 — 9 2 ) > 0. In the case of non-zero surface tension however, 
one only requires a + b + c + d= ^ — r. For this reason (as well as from the pure mathematical 
interest in the analysis of (TlJ) ) , one may as well consider (fijl for all values of the parameters. 



Systems of the form ([T]) have been the subject of intensive investigation over the last decade. 
In particular, the role of the parameters a, b, c, d in the actual fluid models has been explored in 
great detail in the original paper [5] and later in [6]. It was argued that only models in the form 
([!]), for which one has linear and nonlinear well-posedness are physically relevant. We refer the 
reader to these two papers for further discussion and some precise conditions, under which one 
has such well-posedness theorems. 

Regarding explicit traveling wave solutions, Chen, has considered various cases of interest 
in [7], [8]. In fact, she has written down numerous traveling wave solutions (i.e. in the form 
(r/, u) = ((f(x — wt), ip(x — wt)), where in fact some of them are not necessarily homoclinic to zero 
at ±00. In a subsequent paper, [9], Chen has also found new and explicit multi-pulsed traveling 
wave solutions. 

In [11], Chen-Chen- Nguyen consider another relevant case, namely the BBM system, which 
(a = c = 0,b = d= g). They construct periodic traveling wave solutions for the BBM case, as 
well as in more general situations. In [2], the authors explore the existence theory for the the 
BBM system as well as its relations to the single BBM equation. 
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(1) 



Vt + Ux + (rju) x + au x , 

XXX 



- bi] xxt = 
du xxt = 0. 
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(2) 



We wish to discuss another aspect of Q, which is its Hamiltonian formulation. Since it is 
derived from the Euler equation by ignoring the effects of the dissipation, one generally expects 
such systems to exhibit a Hamiltonian structure. This is however not generally the case, unless 
one imposes some further restrictions on the parameters. Indeed, if b = d, one can easily check 
that 

H(r], u) = j —en 2 — au 2 x + rj 2 + (1 + n)u 2 dx 

Furthermore, H (77, u) is positive definite only if a, c < 0. From this point of view, it looks natural 
to consider the case b = d and a, c < 0. In order to focus our discussion, we shall concentrate 
then on this version 

m + u x + (v u )x + au xxx - bn xxt = 
ut + Vx + uu x + crj xxx - bu xxt = 0. 

We will refer to ^ as the Boussinesq 'abc' system. It is a standard practice that stable coherent 
structures, such as traveling pulses etc. are produced as constrained minimizers of the corre- 
sponding (positive definite) Hamiltonians, with respect to a fixed conserved quantity. In fact, 
this program has been mostly carried out, at least in the Hamiltonian cases, in a series of papers 
by Chen, Nguyen and Sun. More precisely, in [12], the authors have shown that traveling waves 
for exist in the regime \b = d, a, c < 0, ac > b 2 . In addition, they have also shown stability of 
such waves in the sense of a 'set stability' of the set of minimizers. In the companion paper |13j . 
the authors have considered the general case b = d > 0, a, c < 0, which in particular allows for 
small surface tension. 

The existence of a traveling wave was proved for every speed \w\ 6 (0,min(l, ^^)).This is the 
so-called subsonic regime. Finally, we point out to a recent work by Chen, Curtis, Deconinck, 
Lee and Nguyen, [10J in which the authors study numerically various aspects of spectral stabil- 
ity/instability of some solitary waves of 0, including the multipulsed solutions exhibited in 
In the same paper, the authors also study (numerically) the transverse stability/instability of the 
same waves, viewed as solutions to the two dimensional problem. 

The purpose of this paper is to study rigorously the spectral stability of some explicit traveling 
waves in the regime b = d > 0, a, c < 0. This would be achieved via the use of the instabilities 
indices counting formulas of Kapitula, Kevrekidis and Sandstede, [15], [16] and the subsequent 
refinement by Kapitula, Stefanov |17j . 

1.2. The traveling wave solutions. In this section, we follow almost verbatim the description 
of some explicit solutions of interest of Q, given by Chen, [7J, see also the more detailed exposition 
of the same results in [8]. More precisely, the solutions of interest are traveling waves, that is in 
the form 

7] = (p(x — wt), u(x, t) = ij)(x — wt). 
A direct computation shows that if we require that the pair (99, ip) vanishes at ±00, then it 
satisfies the system 

(1 + cd 2 )<p - w(l - bd 2 )*P + = 
-io(l - bd 2 )ip + (1 + ad 2 x )^j + ^ = 0. 

The typical ansatz that one starts with, in order to simplify the system Q to a single equation 
is if) = Bip. This has been worked out by Chen, [7], [8]. The following result is contained in the 
said papers. 



(3) 



^which in particular requires that a + b + c + d<0, corresponding to a "large" surface tension r > | 
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Theorem 1. ( Chen, [7], [8] ) Let the parameters a, b, c in the system satisfy one of the following 

(1) a + b £ 0, p = f±| > 0, (p - l/2)((6 - a)p - 6) > 

(2) a = c = -b, b > 

Then, there are the following (pair of) exact traveling wave solutions (i.e. solutions of ([3 
(<p(x — wt), ip(x — wt)), where 

(p(x) = 7]Qsech 2 (Xx) 
ip(x) = B(r)o)r]osech 2 (\x) 

and 



3 + 27/Q 1 2r/ 
w = w{r] ) = ± —====; A = -\ —, , , — r] B{r) 



V3(3 + r/ )' 2y3(a-6) + 26(r/ + 3)' WU; V% + 3' 
and ryo is a constant that satisfies 

(1) r/o = m Case (1) 

(2) 7/o > -3,7/0 / in Case (2). 

1.3. Different notions of stability. Before we state our results, we pause to discuss the various 
definitions of stability. First, one says that the solitary wave solution (<p w , ip w ) is orbitally stable, 
if for every e > 0, there exists 5 > 0, so that whenever ||(/, g) — (<pw,ipw)\\x < ^> one has that 
the corresponding solutions (rj,u) : (/, g) = (rj,u)\t=o 

supinf 11(77(2; — xo,t),u(x — xo,t)) — (tp(x — wt),ip(x — wt))\\x < £■ 
t>o x ° 

Note that we have not quite specified a space X, since this usually depends on the particular 
problem at hand (and mostly on the available conserved quantities), but suffices to say that 
X is usually chosen to be a natural energy space for the problem. This notion of (nonlinear) 
stability has been of course successfully used to treat a great deal of important problems, due 
to the versatility of the classical Benjamin and Grillakis-Shatah-Strauss approaches. However, 
it looks like these methods are not readily applicable (if at all) to the Boussinesq 'abc' system. 
We encourage the interested reader to consult the discussion in [12], where a weaker, but related 
stability was established in the regime ac > b 2 and additional smallness assumption on the wave 
is required as well. This is why, one needs to develop an alternative approach to this important 
problem, which is one of the main goals of this work. 

In this paper, we will concentrate on spectral stability. There is also (the closely related and 
almost equivalent) notion of linear stability, which we also mention below. In order to introduce 
the object of our study, as well as to motivate its relevance, let us perform a linearization of the 
nonlinear system Q. Using the ansatz 

•q = ip{x — wt) + v(t, x — wt) 
u = tp(x — wt) + z(t, x — wt), 

m (|2j) and ignoring all quadratic terms in the form 0(v 2 ),0(vz),0(z 2 ) leads to the following 
linearized problem 

n M 2 ) (v\_ ! + «^ bwd 2 x + ^-w\ 

V- M *>\ z ) t ~ x \ 1 ) \ bwd% + ip-w l + adl + ip J 

Letting 

/ 1 + cd 2 bwd 2 + ij-w\ 2 ! / 1 

(4) L -\bwd 2 + ^-w l + ad 2 + <p 1>J- Ml l*>z) I 1 o 
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the linearized problem that we need to consider may be written in the form 

(5) ut = JLu 

Note that in the whole line context, L is a self-adjoint operator, when considered with the natural 
domain D(L) = H 2 (R}) x H 2 (R}). Letting H := JL, we see that the problem ^ is in the form 
ut = Hu. The study of linear problems in this form is at the basis of the deep theory of Cq 
semigroups. Informally, if the Cauchy problem ut = Hu has global solutions for all smooth 
and decaying data, we say that H generates a Co semigroup {T(t)}t>o via the exponential map 
T(t) = e . Furthermore, we say that we have linear stability for the linearized problem ut = Hu, 
whenever the growth rate of the semigroup is zero or equivalently lim^oo e _5< ||T(t)/|| = for all 
5 > and for all sufficiently smooth and decaying functions /. Finally, we say that the system 
is spectrally stable, if cr(H) C {z : "Rz < 0}. It is well-known that if H generates a Co semigroup, 
then linear stability implies spectral stability, but not vice versa. Nevertheless, the two notions 
are very closely related and in many cases (including the ones under consideration), they are 
indeed equivalent. For the purposes of a formal definition, we proceed as follows 

Definition 1. We say that the problem ^ is unstable, if there is f E H 2 (R}) x H 2 (R}) and 
A : 5RA > 0, so that 

(6) JLf = Af. 

Otherwise, the problem ^ is stable. That is, stability is equivalent to the absence of solutions of 
([§ with A : 5ftA > 0. 

1.4. Main results. We are now ready to state our results. We chose to split them in two cases, 
just as in Theorem [TJ For the case a = c = —b, b > 0, we have 

Theorem 2. Let a = c = —b, b > 0. Then, the traveling wave solutions of the 'abc' system 



2 i x — wt\ . / 3 2 I x ~ w t 



(7) V° SeCk {-2Vf) > ±r >ivoT3 SeCh V 2Vb 

with speed w = ± 3+2r ? Q = are stable, for all r/o '■ r]o 6 ( — 1>0). Equivalently, all waves in (171) are 
a /3(3+t?o) * 1 - r 

stable, for all speeds \w\ < 1. 

Note that \w\ < 1 is equivalent to t]q G (— f,0), so we assume this henceforth. In the remaining 
case, we assume only a = c < 0,b = d > 0, but observe that in this case, Theorem [T] requires that 
r/0 = —3/2, w = 0, that is the traveling waves become standing waves. 

Theorem 3. Let a = c < 0,6 = d > 0. Then, the standing wave solutions of the Boussinesq 
system 

(p(x) = — sech 2 ( — ; — 1 ,tb(x) = ±—=sech 2 
are spectrally stable if and only if 

(8) ((ad 2 x + 1 - p)- 1 ^ - b<p"), (<p - V)> < (I + 12 l ' ' hl 



v 2 ' 5 \a\ 10 a 2 

In particular, the condition Q holds ( and thus the waves are spectrally stable), whenever 

b 

< — < 8.00163, 
—a 
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On the other hand, the condition Q fails ( and thus the waves are spectrally unstable), if 

b 

— > 8.82864. 
—a 

Remark: Note that while, we cannot explicitly compute the value (ad 2 + 1 — (/j) _1 (</? — bip") 
in Q, we obtain estimates, which imply some pretty good results for the stability/instability 
intervals. One can in fact push this further to narrow the gap between the stability and instability 
regions, predicted by Q. This can be done in principle with any degree of accuracy, but it 
increases the complexity the argument. 

2. Preliminaries 

In this section, we collect some preliminary results, which will be useful in the sequel. 

2.1. Some spectral properties of L. We shall need some spectral information about the op- 
erator L. We collect the results in the following 



Proposition 1. Let a, c < and w : < \w\ < min ^1, ^pj • Then, the self-adjoint operator L 
has the following spectral properties 

• Then the operator L has an eigenvalue at zero, with an eigenvector f ^ 

• There is ft > 0, so that the essential spectrum is in a ess {L) C [ft, oo). 

Proof. The first property is easy to establish, this is the usual eigenvalue at zero generated by 
translational invariance. For the proof, all one needs to do is take a spatial derivative in the 

# ) y o 

Regarding the essential spectrum, we reduce matters to the Weyl's theorem (using the vanishing 
of the waves at ±oo), which ensures that 



defining system ([3]), whence L 



o-ess.(L) = o- ess .(Lo) = a\ 



1 + cd 2 bwd 2 — w 
bwd 2 — w 1 + ad' 2 ;. 



That is, it remains to check that the matrix differential operator Lq > k. By Fourier transforming 
Lq, it will suffice to check that the matrix 



1 - c£ 2 -w(b^ + 1 
w(b£ 2 + 1) 1 - < 2 



is positive definite for all £ £ R . Since 1 — c£ 2 > 1, it will suffice to check that the determinant 
has a positive minimum over £ G R 1 . We have 

det(L (f )) = £ 4 (« c - &2 ^ 2 ) + f(~ a - c - 2bw2 ) + C 1 - w 2 ) > (! - w 2 ) + 2 f{Vac - \b\w 2 ), 

where in the last inequality, we have used —a — c > 2^/ac. The strict positivity follows by 
observing that ^fac > \b\w > \b\w 2 , since w < 1. □ 

2.2. Instability index count. In this section, we introduce the instability indices counting 
formulas, which in many cases of interest can in fact be used to determine accurately both 
stability and instability regimes for the waves under consideration. As we have mentioned above, 
this theory has been under development for some time, see |18j . [13] , [19], but we use a recent 
formulation due to Kapitula-Kevrekidis and Sandstede (KKS), [15] (see also [E]). In fact, even 
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the (KKS) index count formula is not directly applicable^] to the problem of ([5]), which is why 
Kapitula and Stefanov, [17] have found an approach, based on the KKS of the theory, which 
covers this situation. In order to simplify the exposition, we will restrict to a corollary of the 
main result in |17j . More precisely, a the stability problem in the form is considered in the form 

(9) d x Cu = An, 

where £ is a self-adjoint linear differential operator with domain D(C) = H S (R ) for some s. It 
is assumed that for the operator C, 

(1) there are n(C) = N < +oo negative eigenvalues^] (counting multiplicity), so that each of 
the corresponding eigenvectors {fj}jLi belong to ff 1//2 (R 1 ). 

(2) there is a k > such that a ess (C) C [k 2 , +oo) 

(3) dim[ker(£)] = 1, ker(£) = span{ipo}, ipo real-valued function, Vo £ H^CR 1 ) nif _1 (R 1 ). 
Here, if _1 (R 1 ) is the homogeneous Sobolev space, defined via the norm 

: = (/ R1 ^p 1 ^) • 

or equivalently, u = d x z in sense of distributions, where z£L 2 and := IMIl 2 - I n that 

case, we have 

Theorem 4. (Theorem 3.5, [T7] ) For the eigenvalue problem 

(10) d x Cu = Xu, nGL 2 (R 1 ), 

where the self-adjoint operator £ satisfies D(C) = ii" s (R ) for some s > 0, assume that 

Then, the number of solutions of n unsta bi e (£) , with A : 3?A > safe/ie^] 

(11) < n unsta ble(d x £) = n(C) - n ({C~ x d~ x ^ Q , d~ Vo)) mod 2. 

Of course, our eigenvalue problem ^ does not immediately fit the form of Theorem [4] First, 
Theorem [4] applies for scalar- valued operators £, while we need to deal with vector- valued opera- 
tors. This is a minor issue and in fact, one sees easily that the arguments in |17| carry over easily 
in the case, where £ is a vector-valued self-adjoint operator as well. A second, more substantive 
issue is that the form of Q is not quite the one in (10). Namely, we have that the operator J, 
while still skew-symmetric is not equal to d x . 

In order to fix that, we need to recast the eigenvalue problem ^ in a slightly different form. 
Indeed, letting f = (1 — 6c3 2 ) _1 / 2 g and taking (1 — bd x ) 1 / 2 on both sides of we may rewrite 
it as follows 

l ) a - Mgr 1/a £(i - ^ 2 r 1/2 g = Ag. 

J J); L := (1 - bdl)-^L{l - 6fl£)-V> 



If we now introduce 



2 due to a crucial assumption for invertibility of the skew-symmetric operator J, which is not satisfied for d x 
acting on R 1 

^We will henceforth denote by n(M) the number of negative eigenvalues (counting multiplicities) of a self-adjoint 
operator M 

^here d x 1 %[)o is any L 2 function /, so that = d x f in distributional sense 
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we easily see that J is still anti-symmetric, L is self-adjoint and we have managed to represent 
the eigenvalue problem in the form JLg = Ag. Note that the operator J is very similar to d x , 

except for the action of the invertible symmetric operator f ^ J ~\ on it. It is not hard to see 

that the result of Theorem [4] applies to it (while it still fails the standard conditions of the KKS 
theory, due to the non-invertibility of J). Note that one needs to replace d~ l by J -1 in the 
formula (11). Furthermore, the number of unstable modes for the two systems {JL and JL) is 
clearly the same, due to the simple transformation (1 — bd 2 )" 1 / 2 connecting the corresponding 
eigenfunctions. 

Thus, if we can verify the conditions under which Theorem^ applies, we get the stability index 
formula 



L ( ^ J = 0, we conclude that L[(l - bd 2 ) 1 / 2 ( ^ )] = 0. It follows 



(12) n unstab ie(J L) = n unstable (J L) = n(L) - n({L 1 J 1 i/j ,J Vo)) mod 2. 

Since by Proposition 
that Vo = d x (l - bd 2 ) 1 / 2 J and 

(Z^j-W-Vo) = (L-'Kl - bd 2 x ) (I J)(j)],(l-Mg 



lUMwi ,W 1 W <p 

io U 



<irM(i-^)(£)],(i-w£)(£): 



Thus, we conclude that we will have established spectral stability for ([6]), if we can verify the 
conditions (1), (2), (3) of Theorem kl for the operator L, n(L) = 1 and 



(13) (L-^l-bdl)^ )},(l-bd 2 x )^ ))<0. 

and instability otherwise. 

Concretely, we will verify the conditions on L in Proposition [2] below, after which, we compute 
the quantity in (13) in Proposition [3j 

Proposition 2. The self-adjoint operator L = (1 — bd 2 )~ l l 2 L(l — bd 2 )~ 1 / 2 satisfies 

(1) a ess ,(L) C [k, 00) for some positive k. 

(2) n{L) = 1. 

(3) Ker(L) = span{(l - Mg)V» ( ^ )}. 
in ifoe following cases 



.a = c = -b,b>0,B = ±^ q> W = ±-^M= y % e (-1,0). 
• a = c<0 ; 6>0,u> = 0,5 = ±>/2- 
Proposition 3. Regarding the instability index, we have 

. For a = c = -b, b > 0, w = ± Jgg^ , B( Vo ) = ±J^I, and for all m € (-|,0), 

-).r, , / V» n /1 ^ f ^ 



<L- A [(l-6^)( J )],(l-&fl£)( J ))<0 
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For a = c < 0, b > 0, w = 0, B = ±y/2, 



/„, — / 9 12 6 3 2 

- — 

3 V \ 2 5 a 



in particular, 



(L- l [(l-bd 2 x )[ t )],(l-bd%)( t )><0, 0< — < 8.00163, 



<P J \ <P J -a 

{L- l [(l (1 - 6fl£) ))> 0, — a > 8.82864. 

Theorem [2] follows by virtue of Proposition [2] and Proposition [3} Thus, it remains to prove 
these two. 

3. Proof of Proposition [2] 
We start with the gap condition for a ess ,(L) stated in Proposition [2] 

3.1. L is strictly positive. The idea is contained in Proposition [l] Write 

l = {i-bdir^L{i-bdir^ = (i-65 x 2 )- i / 2Lo(1 _^ ) -i/ 2+( i_^ r i/ 2(L _ Lo)(1 _ 6a 2 r i/ 2i 

where L — Lq is a multiplication by smooth and decaying potential. It is also not hard to see 
that (1 — bd 2 )^ 1 ^ 2 is given by a convolution kernel K : K(x) = f°° , 6 = =d£, which decays 

faster than polynomial at ±00. It follows that the operator (1 — 6(9 2 ) _1 / 2 (L — £o)(l — is 
a compact operator on L 2 (R X ) and hence By Weyl's theorem 

o ess .(L) = a ess .((l - &<9 2 r 1/2 Lo(l - ftflg)- 1 / 2 ) = tr((l - 69 2 )- 1/2 A>(l - bS 2 )- 1 ^) 
Thus, as we have explained in the proof of Proposition [TJ it will suffice to check that the matrix 

(1 + 4vr 2 6c 2 )- i / 2 L (o(i + ^ner 1 ' 2 

is positive definite. But since Lq{£) is positive definite, the result follows. Note that this only 
shows that a ess ,{L) > 0. Since we need to show an actual gap between a ess ,{L) and zero, it suffices 
to observe (by the arguments in Proposition [l]) that the eigenvalues of Lq(£) have the rate of 0(£ 2 ) 
for large £, which implies that the positive eigenvalues of (1 + 47r 2 &£ 2 ) -1 1 2 Xo(£)(l + 4-7r 2 6£ 2 )~ 1 / 2 
have the rate of O(l). 

3.2. The negative eigenvalue and the zero eigenvalue are both simple. We now pass to 
the harder task of establishing the existence and simplicity of a negative eigenvalue for L as well 

as the simplicity of the zero eigenvalue. Note that as we have already observed L f "\ =0. It 

follows that 

m - bdlf' 2 ($)] = (!- bd 2 x )-^ 2 [L ( £ )] = 0. 

Thus, we have already identified one element of Ker(L), but it still remains to prove that 
dim{Ker{L)) = 1, in addition to the existence and the simplicity of the negative eigenvalue 
of L. 
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Next, we find it convenient to introduce the following notation for the eigenvalues of a self- 
adjoint operator £. Indeed, assume that £ = £* is bounded from below, £ > — c, we ordeij^the 
eigenvalues as follows 

inf spec(£) = Xq(£) < Ai (£)<... . 
Recall also the following max min principle, due to Courant 

A (£) = inf (£/,/), Ax (£) = sup inf (Cf,f), A 2 (£) = sup inf <£/,/>. 

Il/ll =1 g^O ll/ll=lj/-Lfl gi,gr-9l^a92 \\f\\= 1 >f- Ls P an l9l,92] 

Clearly, our claims can be recast in the more compact form 
(14) A (L) <0 = Ai(Z) < A 2 (L). 

matters from L to standard second order differential operators, like L. 

Lemma 1. Let a, c < 0, b > and w : < \w\ < min ^1, ^ygj^ • Then 

• all eigenvectors of L from Q, corresponding to non-positive eigenvalues, belong to 
H°°(R l ) = nfZ 1 H l (R 1 ). 

• If C is any bounded from below self-adjoint operator, for which 

Ao(jC) < = Ai(£) < \2(C), and S is a bounded invertible operator, then 

\ (S*£S) < = X^STCS) < X 2 {S*£S). 

• If L has the property Xo(L) < = Ai(L) < ^(L), then so does 
L = (1 - bd^)-^ 2 L{l - bdl)- 1 ' 2 . That is, ph ZioWs. 



Proof. (Lemma [T]) 

Take the eigenvector f, corresponding to —a 2 , a > 0, i.e. Lf = — a 2 f. As observed in the 
proof of Proposition [TJ we can represent L = Lq + V, where V is smooth and decaying matrix 
potential. In addition, recall Lq > k, hence Lq + a 2 > nld and hence invertible. It follows that 
the eigenvalue problem at —a 2 can be rewritten in the equivalent form 

f = -(L + a 2)-i[ V f] 

Clearly, (Lo+a 2 ) -1 : L 2 — > H 2 , whence we get immediately that / G H 2 , if / 6 L 2 . Bootstrapping 
this argument (recall V G C°°) yields / G H 4 , H 6 etc. In the end, f G H°°. 
Next, we have 

MS'CS) = ^UrCSf, f) = inf = inf ^ < 0, 

since \q(C) = mi g -\\ g \\ = i{Cg , g) < 0. Since Ai(£) = 0, it follows that there is h, so that 
iai g ± h (£g,g) > 0. Thus, 

\ /cvc\s ■ r ( S *£- S fJ) ■ f (£9,9) ^ n 

Xx(S CS) > mf r-—. = mf — > 0. 

f±Sh ll/H 2 g±h WS^gW 2 

Since is still an eigenvalue for C with say eigenvector x, it follows that S~ 1 x is an eigenvector 
to S*CS, so is also an eigenvalue for S*CS and hence X\(S*CS) = 0. 

Regarding X2(S*£S), we already know that X2{S*CS) > Xi(S*£S) = 0. Assuming the contrary 
would mean that X2(S*£S) = 0, that is is a double eigenvalue for S*£S, say with linearly 
independent eigenvectors fi, f 2 - From this and the invertibility of S, it follows that S~ l f\, S^ 1 
are two linearly independent vectors in Ker(L), a contradiction with the assumption that is a 
simple eigenvalue for L. 

^We follow the standard convention that if an equality appears multiple times in the sequence of eigenvalues, 
that signifies that eigenvalue has the same multiplicity 
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The result regarding (1— bd 2 )~ l l 2 L(l — bd 2 )^ 1 / 2 follows in a similar way, although clearly cannot 
go through the previous claim (since (1 — bd 2 )^ 1 ^ 2 does not have a bounded inverse). To show 
that Xq{L) < 0, take an eigenvector say go : \\go\\ = 1, corresponding to the negative eigenvalue 
—a 2 for L. Note that by the first claim, such a go is smooth, so in particular (1 — bd 2 ) l l 2 go is 
well-defined, smooth and non-zero. We have 



(L(l - bd^go, (1 - bd^go) _ (Lg ,g 



2 



° { ' ~ 11(1 " bd*)V*go\\* ~ ||(1 - bd 2 y/ 2 g \\ 2 ~ 11(1 - bd*) 1/2 9o\\ 2 ' 

Next, to show that Ai(£) > (the fact that is an eigenvalue for L was established already), 
recall that since L has a simple negative eigenvalue, with eigenfunction go, we have 



inf (Lg,g) = 0. 

9:9-1-30 



It follows that 



wn^ -f (Mill - ■ f W*) 

Ai(L) - /k4-^ ttf - i nd - 

Regarding the proof of \2{L) > 0, we start with \2(L) > Xi(L) = and we reach a contradiction 
as before (i.e. we generate two linearly independent vectors in Ker(L)), if we assume that 
A 2 (Z) = 0. □ 

Using Lemma [TJ allows us to reduce the proof of ( 14 ) to the proof of 

(15) Ai(L) <0 = Ai(L) < A 2 (L), 

which we now concentrate on. 
We have 



1 + ad 2 bwd 2 +ip — w 
bwd 2 + ijj — w 1 + ad 2 + <p 

{l + adl)Id+{bwdl- W )( { \ + ( ° ' 



L 



1 J \ ip <p 

( — —\ 

Introduce an orthogonal matrix T = ( *\ J and observe that 



V2 V2 

I u J = T l -1 



° ^ T-Ul ° IT. 



It follows that 



L = T- 1 ((l + ad 2 )Id+(b W d 2 - W + ^( I [jjl, 

whence, by unitary equivalence, it suffices to consider the operator inside the parentheses. That 
is, we consider 

nm M _( -dl(-a-bw) + {l-w)+^ + % % 
[Lb) M ~\ f -flg(-o + 6«;) + (l + «j)-^ + | 

We shall need the following 
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Lemma 2. Let a, A > and Q £ R . Then, the Hill operator 

£ = -d 2 x + a 2 - Qsech 2 (Xx) > 

if and only if 

(17) a 2 + aX>Q. 

Proof. This is standard result, which follows from the ones found in the literature by a simple 
change of variables. First, if Q < 0, we see right away that C > and also the inequality (17) is 
satisfied as well. So, assume Q > 0. Consider Cf = of and introduce f(x) = g(Xx). We have 
(after dividing by A 2 and assigning y = Xx) 

r Ct 2 Q O/ M ° I \ 

l-°yy + T2 - J2 sech W)\9 = ^{y) 
Recall that the negative the operator —d yy — Zsech 2 (y) are k m = 



m 



, pro- 



vided (Z + I) 



m 



2 > 0, m = 0, 1, 2... [see pQ]. Note that ko = inf o~(—d yy — Zsech 2 (y)) and 



hence, to avoid negative spectrum, we need to have 



a 



0<^ + k 



of 
A 2 



9l , I v 

A 2 4 



Solving this last inequality yields (17). 



□ 



We are now ready to proceed with the count of n(L) in each particular case of consideration. 
Case I: a = c = —b, b > 

Going back to the operator M, we can rewrite it as 



where S 
to 




_£_ 







. Thus, according to Lemma 



we have reduced matters 



1. 



M 1 = (-d 2 + -)Id + <p 



2b^l-w 2 
b(l+w) 

Diagonalizing this last symmetric matrix yields the representation 



b(l-w) 
I 

2b\/\-w 2 



B-\ 



b(l-w) 2bVl-w 2 



( l+2Bw+V4B 2 +4Bw+l 
2b(l-w 2 ) 
q l+2Bw-V4B 2 +4Bw+i 





2b(l-w 2 ) 



u 



for some orthogonal matrix U. Factoring out U* , U again and using Lemma [T] once more reduces 
us to the operator 



M 2 



Ci 

c 2 
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which contains the following Hill operators on the main diagonal 

o2 , 1 1 + 2Bw + V4B 2 + 4Bw + 1 u2 

L\ = -o x + - + r/o —r-y- Y\ sec/l 



b ,u 26(1-^2) \ 2 y/b 

X 



, 1 1 + 2Bw - yjAB 2 + 4Bw + 1 , 9 



+ I + ^0 T^Tn sech 

Using the formulas 



b ,u 26(i-^2) y 2y /b 

Note that n{L) = n(d) + n(C 2 ). 



yields 



1 3 

Ci = -dl + - — -sech 2 , _ , . 

b & V2vV 

£ 2 = 3*70 , ec/t 2 



6 6(9 + 4%) \2Vb 
According to the formulas for the eigenvalues in Lemma |2j (with a = A = i^j^iQ = f > 0) we 
have that ^ 

MA) = ^ - ( l/f^| - §) =2- ( v / 1^25-1.5) 2 = 0, 



which indicates that £i has one negative eigenvalue and the next one is zero, whence n(C\) = 1 

9 

4' 



for all r]o > —3. Thus, n(L) = 1 + n(C_)- It is also immediately clear that for 770 £ (— 1,0), 



£2 > and hence n(L) = 1. 

Case II: a = c<0,6 = d>0, a + b^0 

In this case, we have p = = 1, 7/0 = 3 ^ 1 2p 2p ^ = — | and thus w(t]q) = w(—3/2) = 0, 
A = 2v e^ ' B(rjo) = ±y/2. This simplifies the computations quite a bit. In fact, starting from the 
operator M, defined in (16), we see that it has the form 

M = {adl + l)Id+(^ B f^ ^ 

Recall that here B = ±y2- Consider first B = \[2. Diagonalizing the matrix vian an orthogonal 
matrix S yields the representation 

V2 + 1 \ \ _ _( 2 \ 

J_ / ^3 + 2^/2 ^3-2^2 \ 
" \/6 I -V3- 2v/2 ^3 + 2^ J 



S 



Thus, in this case, we have represented the operator L in the form 

(18 > L ^s T) .(-<^^ _j.lt _ v ) sr. 

where S, T are explicit orthogonal matrices. It is now clear that since r/o = — § < 0, we have that 
<p(x) < and hence the operator ad 2 + 1 — tp > 0. On the other hand, L^dv = + 1 + 2y? 
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is well known to have a zero eigenvalue (with eigenfunction <//) and an unique simple negative 
eigenvalue. 

For the case B = —y/2, we have (18), with 



S 



J_ / \/3-2V2 ^3 + 2^/2 
V§ I a/3 + 2V2 -a/3 - 2 a/2 



4. Proof of Proposition [3] 



The purpose of this section is to compute the quantity appearing in ( 13 ), whose negativity will 
be equivalent to the stability of the waves. Thus, we need to find 



L-\{l-bdl) 



Here, our considerations need to be split in two cases: a = c = —b, and o = c < 0, b > 0. 

The case a = c = — b is easier to manage, since in int we have a a free parameter w = w(rjo) that 
we can differentiate with respect to in ([3]). The remaining case is harder, because the parameter 
770 = —3/2, whence w = and one cannot apply the same technique. 

4.1. The case a = c = —b, b > 0. Taking a derivative with respect to w in ([3]), we find 



d w ip 



whence 



We obtain 



d w ip 
L-i[(l-bdl 



(1 - bd 2 x ) 



d w ip 



{{l-bdl 



d w [(<p,^) + b{<p',ip')} = d w [B( m ) j ^(£) 2 + %W^] = 

Bd w [[ [ V 2 (o+bv ,2 mdi]+d w B I ^ 2 (e)+v 2 (£M 

JR JR 



16Vb 
5 



\<o^fb 
5 



2S + 770 



dr] 



dr] , . 

m-r- =: d ( w ) 

aw 



We are now ready to compute this last expression in the cases of interest. 



d w (p 
d w ip 



4.1.1. B(rj ) = -J^- t w = , 3+2r?0 . We have 



d??o 2^/3(3 + 7/0)1 dB a/3 1 



3 



2?7o + 9 ^770 2 ( 3 + ??0 ) 2 
and 



48 V 36 N drjo 

ayai) = ' 

1C 

for -| < r/o < 0. 



d(w) = r (4 + Vo)vo-f- < 

10(3 + 7/o)5 aw 
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4.1.2. 5(770) = Jxg-,w= f +2m . We have 



drjo 2^3(3 + r7o)i ^3 1 



dw 2r/ + 9 ' dr/o 2 (3 + 770)! 

hence 



for -f < r? < 0. 



. 48v36 . drjo 

d(w = — -3- (4 + 770)770 < 0. 

10(3 + 770)2 dw 



4.2. The case: a = c < 0, b > 0. As we have discussed above, we have explicit formulas for all 
the quantities involved. Namely, we have w = 0, A = i^= a ,B = ±y/2. Thus, 



<p(x) = — -sech 2 



2 V2^/=^ 

4.2.1. Case B = y/2. We need to compute 

( (l-bd 2 x )ijj \ f ( I -M'i)r 



1 - bd^ J ' V (1 - bdl) V 



To that end, we use the representation (18). We have 
x / (1-63^ \ / (1-^ 



H ( °« + X + 2 * adl + \ _ „ ) [sr ( f ) (i - sr ( f ) (i - 

A direct computation shows that ST ^ ^ ^ = ^ ^ , whence our index / can be computed 

as follows 

I = \{{ad 2 x + 1 + 2v)- l [(l - b%)<p],(l - b%)<p) + \{{adl + 1 - <p)-\(\ - b%)<p], (1 - bd 2 x )tp) 
Denote / = (1 — bd^)ip and 

LRdV = ad 2 + 1 + 2y? 
Lhui = ad 2 + l-(p 

Note that by Weyl's theorem a ess .(LHUi) = [l 5 oo). On the other hand, by the fact that tp < 0, 
the potential —92 > and hence, by the results for absence of embedded eigenvalues, a (Lhui) = 
cr ess ,(LHUi) = [l>oo). We now compute the index 

i = \m K i dV f,f) + (L H \ u fj)). 

To that end, we differentiate the equation 

a<p" + ip + ip 2 = 

with respect to a. We gelQ 

(19) L K dV<fa = -¥>", 



°we use the notation y> a = d a ip denotes the derivative with respect to a 
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whence Lj} dV [p"] = —p a . Using that LxdV'P = p 2 = —ap" — ip and the above relation, we obtain 
that 

-<p = aL^yip" + L-^ytp = -cupa + L^yif. 

It follows that 

(20) ^KdvV = a( Pa~P, 

(21) L~ K \yf = ( a + b)p a -p. 
and 

( L /fdv/' /) = (« + b )(p a , p) ~ b(a + b)(<p a , p") - (p, p) + b(ip, if/'). 
By direct computations 



1 d [ + °° 2 3 
- — / <p ax = — , 

2 da /_ no 2J—a 



/+oo 
-oo 



{p a ,p") = ~\~T~ I p' 2 dx = - : 



2daJ_ nci lOlalv^a 



(<£>,<£>) = ^V^a I sech 4 (y)dy = 6^J^a. 



As a consequence, 



X K *v J,J/ 2^a 10\a\y/=a h^a 



9 , — 12 6 36 2 , — / 9 12 b 3 b 2 



2 5 ^ lOlalv 73 ^ V 2 5 |o| 10 a 2 , 

This yields the desired computation for the terms involving Lj} dV . We turn our attention 
to Ljj 1 ^. The situation here is a bit trickier, since we cannot compute explicitly the quantities 
Lff\n [<p] , [ip ff ] , as required in the formula for /. Instead, we need to rely on estimates. To 
start with, observe that 

Lmu[p] = ap" + p- p 2 = -2p 2 = 2ap>" + 2p, 

whence 

(22) + ¥>] = !■ 

Since we need to compute LJj\ u [f] = L~^ iU [p — bp"], we project the vector / onto ap" + p and 
its orthogonal subspace as follows 

(p - bp", ap" + p) 
\ap" + p\ 

Calculations then show that since 



f = p-bp = , , „ 2 {ap +p) + g 



I //i|2 / // Ih ( '' 

\<P II = \f ,<P 



7\a\y/—a 
we have that 



'=(Kr) (V ' +v)+9 
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whence 

£ ™W-(5 + §jHj)5 + £ ™M- 

Thus, the quantity that needs to be computed is 



1/7 2 6 
2 U 



2 6 \ / . // V ,r-l V 1 2 6 A 

+ 9hJ ^-^^ ) + (L ^ 5) + 2U + 9hJ 

All of these can be computed explicitly, except for (L~^ iU g , g) , which we estimate by 
< (L~^ iU g,g) < \\g\\ 2 , which holds since (t{Lhui) C [1,oo). Thus, 



22 6 2 2 6 26 
45 + 9 y + "9" 



and on the other hand 



1 /7 2 6 \ . . , — / 4 6 2 46 6 112 

2 9 + 9R ( ^ ) = ^te + «R + ¥ 



Thus, we obtain the following estimate for the instability index I 

- - S ( ^,/) + (^/,/)^(s(-^-^ + £) + (|^^ + f 



2a/^ / 6 2 6 

65 — — 427 — 1= — 745 



45 \ a 



,2 



'—a 



On the other hand, we have the following estimate from below 

3, - s^U/) + (^/./)>^(s(-^-^ + ] £) + (^ + ^ + f ) 



6 2 6 
.56^-409^- 754). 
45 \ a 2 \a\ 



The picture below shows the graphs of the two estimates of 31/y/^a. If one solves the corre- 
sponding quadratic equations, we see that we have stability, whenever 

< — < — ( 427 + 3741781) ~ 8.00163. 
—a 130 V / 

and instability, when 

— > — (a09 + 3737353^ ~ 8.82864. 
—a 112 V / 
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20 




Figure 1. The picture shows the graphs of the function ^ (65z 2 - 427z - 745) 
is in blue, while ^(56z 2 — 409z — 754) in red. Note that the graphs do coincide 
for z = 1, which is the case, since there g = and the computations becomes 
precise. 

4.2.2. Case B = —y/2. In this case, the computation for the index is the same since 
r (T -i( \ / (l-bd^ \ 

= <(^™ ^ + \_ v )^(-f )d-^],5T(-f )(!-«©*,> = 
= ^(8(^y/>/) + (^/ ! /))> 
where in the last line, we have used that ST ( \ = ( V/^ ^ as above. The rest of the 



1 J L , 

argument proceeds in exactly the same way, since the exact same quantity is being computed. 
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